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ABSTRACT 

We investigate the evolution of the metallicity of the intergalactic medium (IGM) 
with particular emphasis on its spatial distribution. We propose that metal enrichment 
occurs as a two step process. First, supernova (SN) explosions eject metals into 
relatively small regions confined to the surroundings of star- forming galaxies. From a 
comprehensive treatment of blowout we show that SNae by themselves fail by more 
than one order of magnitude to distribute the products of stellar nucleosynthesis over 
volumes large enough to pollute the whole IGM to the metallicity levels observed. 
Thus, a additional (but as yet unknown) physical mechanism must be invoked to mix 
the metals on scales comparable to the mean distance between the galaxies which are 
most efficient pollutants. From this simple hypothesis we derive a number of testable 
predictions for the evolution of the IGM metallicity. Specifically, we find that: (i) 
the fraction of metals ejected over the star formation history of the universe is about 
50% at z = 0; that is, approximately half of the metals today are found in the IGM; 
(a) if the ejected metals were homogeneously mixed with the baryons in the universe, 
the average IGM metallicity would be (Z) = Jl^/il;, ~ 1/25^© at z = 3. However, 
due to spatial inhomogeneities, the mean of the distribution of metallicities in the 
diffusive zones has a wide (more than 2 orders of magnitude) spread around this value; 
(in) if metals become more uniformly distributed at z ^ 1, as assumed, at z = 
the metallicity of the IGM is narrowly confined within the range Z ~ 0.1 it O.OSZq. 
Finally, we point out that our results can account for the observed metal content of 
the intracluster medium. 



Subject headings: Cosmology: theory - intergalactic medium - quasars: absorption 
lines 



1. Introduction 

Primordial nucleosynthesis enriched gas in the universe with the light elements He, D, and 
Li. It is only when the first galaxies and their stars appeared that heavier elements could be 
synthesized and, in some cases, ejected into the intergalactic medium. In currently popular 
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models of galaxy formation based on hierarchical clustering the first galaxies to form were low 
mass systems with such shallow potential wells that a few supernovae could deposit sufficient 
kinetic energy to expel the entire interstellar medium of the galaxy (Ferrara 1998). In this way, 
these initial episodes of star formation in the universe (sometimes referred to as Population III) 
may have enriched the intergalactic medium (IGM) to an average metallicity Ziqm ~ 10~'^Zq 
(Miralda-Escude k Rees 1997, Nath k Trentham 1997, Gnedin k Ostriker 1997, Ciardi et 
al. 2000b), comparable to that of the most metal-poor stars in the halo of our Galaxy (Ryan, 
Norris, k Beers 1996). 

Most of the metals in the universe were presumably produced in larger collapsed structures 
at redshifts z ^ 10. In general such galaxies were better able to retain the products of stellar 
nucleosynthesis and could therefore achieve the abundance levels observed today. It is also likely 
that some metal-enriched gas escaped into the IGM, but it is far from clear by which process 
and to what degree the metals became distributed over large volumes, far from their production 
sites. These are the questions we consider in the present paper. Recent numerical hydrodynamic 
simulations of large-scale structure formation (Hellsten et al. 1997; Ranch et al. 1997; Zhang et 
al. 1998) have shown that baryons in the universe are distributed in a network structure (the 
"cosmic web") with galaxies located at the high overdensity peaks. This implies that most of the 
volume is occupied by voids filled with gas at, or below, the mean cosmic density; such gas can be 
viewed as the long-searched-for true intergalactic medium. Since the voids are very large, with 
typical dimensions of several Mpc, their pollution by heavy elements produced in SN explosions is 
far from a trivial problem. In any case there must have been a protracted epoch when the metal 
content of the IGM was highly patchy. 

Observationally, our best view of the IGM is still provided by the Lya forest in the spectra 
of distant QSOs (Sargent et al. 1980) and the detection of metals in the forest ranks as one of 
the most significant discoveries made possible by the Keck telescopes (Cowie et al. 1995; Tytler et 
al. 1995). From an analysis of these data Hellsten et al. (1997) and Ranch, Haehnelt, k Steinmetz 
(1997) concluded that the measured column density ratios A'^(C IV)/A^(H I) imply that typically 
[C/H] ~ —2.5 at 2 ~ 3, with a one order of magnitude dispersion in the metallicity of different 
clouds about this mean value.Q At lower redshifts {z = 0.3 — 0.8) Barlow k Tytler (1998) using 
HST/FOS spectra, with reasonable assumptions concerning the ionization correction and the line 
clustering properties, conclude that metallicities are as high as [C/H] ^ — 1.3, roughly an order of 
magnitude larger than the value at z = 2.5. These measurements, however, still refer to overdense 
regions of the universe, traced by Lya clouds with column densities in excess of log A^(H I) = 
14.5 . The situation is far less clear-cut when we turn to the voids — or log A^(H I)< 14.0 — where 
the observations are very challenging even with a 10-m telescope. Two studies have addressed this 
problem, with confiicting results. Lu et al. (1998) applied a stacking technique to nearly 300 C IV 
regions in QSO spectra but still found no composite signal. They interpreted this non-detection 



In the usual notation, [C/H] = log (C/H) - log (C/H)©. 
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as evidence for a highly non-uniform degree of metal enrichment at z ~ 3, with the voids having 
[C/H] ^ — 3.5 . On the other hand, Cowic & Songaila (1998) used a pixel-to-pixel optical depth 
technique to conclude that the average C IV/H I ratio remains essentially constant over the full 
range of neutral hydrogen column densities tested, down to log A''(H I) ~ 13.5 . If this is indeed 
the case, the ejection and transport of metals away from galaxies must have been much more 
efficient than envisaged. Most recently, Ellison et al. (1999, 2000) have re-examined the problem 
and found that both approaches suffer from limitations which had not been properly taken into 
account in previous analyses. In their view, whether there are metals in the voids remains an open 
question. 

On the theoretical side too, only a limited amount of work has been done on this topic. The 
common assumption that supernova driven winds may be able to distribute metals over large 
distances has been shown to be too simplistic and does not stand up to close quantitative scrutiny 
(MacLow & Ferrara 1999; Murakami & Babul 1999; D'Ercole & Brighenti 1999). These studies 
have come to the conclusion that efficient blowout is likely to be inhibited by the galaxy ISM and, 
at least at high redshift where densities are higher, by the pressure of the surrounding intergalactic 
gas (Babul & Rees 1992; Ciardi & Ferrara 1997). Gncdin & Ostriker (1998) numerically simulated 
the IGM enrichment in a CDM+A cosmological model. They found that the mean metallicity 
of the universe at z = 4 is about 1/200 Zq but the variations around this value are large and 
dependent on the overdensity. At this epoch the enrichment was incomplete, with some regions 
of the universe still of pristine composition. Their simulations, however, end at z = 4 so that 
a direct connection with the nearby universe is difficult. In addition, quantitative conclusions 
on the transport mechanism are affected by uncertainties arising from insufficient numerical 
resolution so that, for example, they are unable to resolve the snowplough phase of explosions. 
To alleviate the problem of inefficient blowout Gnedin (1998) developed the idea that metals 
are predominantly transported following merger events, a mechanism which would also predict a 
highly inhomogeneous distribution of heavy elements in the IGM at z = 4. However, the merger 
history has not yet been followed to more recent epochs. 

In this paper we make explicit predictions for the IGM metallicity evolution by suggesting 
a two-step mechanism for the transport and mixing of heavy elements. In our scenario, SN 
explosions first eject metals over a relatively small region in the surroundings of the host galaxy; 
subsequently, a yet unknown diffusion process transports and mixes these elements over much 
larger scales. The plan of the paper is as follows. In §2 we consider the conditions for gas ejection 
from galaxies, paying particular attention to the treatment of blowout, which we show to be quite 
a rare occurrence in massive galaxies. In §3 we explicitly derive the metallicity evolution of the 
IGM for a particular (Cold Dark Matter) cosmological model. The results are discussed in §4 and 
a brief summary (§5) concludes the paper. 
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2. Metal Ejection by Galaxies 

It is likely that the first episodes of star formation in the universe had a dramatic impact on 
the small mass galaxies in which they occurred. The combined effects of supernova explosions 
in relatively small volumes lead to the formation of superbubbles (SBs) which can result in 
partial (blowout) or even complete (blowaway) removal of the interstellar medium (Kovalenko 
& Shchekinov 1985; MacLow & McCray 1988; Ciardi k Ferrara 1997; Ferrara 1998; MacLow 
& Ferrara 1999). The resulting large scale outflows (galactic superwinds) have been observed 
directly in local starbursts (Gonzalez Delgado et al. 1998; Heckman et al. 1998) and in Lyman 
break galaxies at z ~ 3 (Lowenthal et al. 1997; Pranx et al. 1997; Pettini et al. 1998; Pettini et 
al. 2000). For the purpose of this paper, we are interested in establishing the maximum galactic 
mass which allows blowout to occur. MacLow &: Ferrara (1999) have shown that when blowout 
does take place the escape efficiency of the metals produced by the supernova progenitors is close 
to unity. 



2.1. Conditions for Blowout 

The blowout condition can be derived by comparing two characteristic velocities: the blowout 
velocity, Vh, and the escape velocity of the galaxy, Vg. To calculate these two velocities we start by 
defining a protogalaxy as a two component system consisting of a dark matter halo and a gaseous 
disk. We assume a modified isothermal halo density profile Ph{r) = Pc/[^ + (f/i^a)'^] extending out 
to a radius r2oo ^ = [3M/j/47r(200/9crit)]^^^) defined as the radius within which the mean dark 
matter density is 200 times the critical density pcnt = ^Hq{1 + 2;)^/87rG at the redshift z when the 
halo is identified. is the halo mass and Hq = lOO/i km s~^ Mpc~^ is the present-day Hubble 
constant; throughout this paper we adopt = 1. 

For such a halo the escape velocity can be written as 

vl ~ A-KpGpcvl = (1) 

with p = 1.65 (MacLow & Ferrara 1999). Note that Ve = \/2p Vc, where Vc is the circular velocity 
of the halo. 

Due to the dissipative nature of the gas, the baryons which initially should be distributed 
approximately as the dark matter, will lose pressure and collapse in the gravitational field of the 
dark matter. The density ratio between these two components is assumed here to be equal to its 
cosmological value, and therefore the initial density of the gas in the protogalaxy is pg = ^bPh- If 
the halo is rotating, the gas will collapse in a centrifugally supported disk. The radius of the disk 
can be estimated by imposing that the specific angular momentum of the disk, j^, is equal to that 
of the halo, jh (Mo et al. 1998, Weil et al. 1998) 



Jh 



V2Xvcrh = 2vcid = jd (2) 
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where id is the disk scale length and A the standard halo spin parameter. As shown by numerical 
simulations (Barnes & Efstathiou 1987; Steinmetz &: Bartelmann 1995) A depends very weakly 
on M/i and on the density fluctuation spectrum; its distribution is approximately log-normal and 
peaks around A = 0.04. From eq. ^ we then obtain = {\/y/2)rh- For an exponential disk, 
implicitly assumed in deriving eq. |2| above, the optical radius (i.e. the radius encompassing 83% 
of the total integrated light) is 3.2 x H^- Since galaxies typically extend ~ 2 times their optical 
radius in H I (Salpeter &: Hoffman 1996), we adopt for the radius of the gaseous disk 

Td = 4.5Ar;,. (3) 

We now evaluate the scale height of the gas in the disk. This is roughly given by 

where we have used the expression for above; Cg is the effective gas sound speed which also 
includes a possible turbulent contribution. Note that 

- = 15.3A('^V (5) 

rd \VeJ 

If all baryons bound to the dark matter halo were able to collapse, the mean gas density in the 
disk would be 

rh\^ (rd 



The previous assumption, although uncertain, gives approximately the correct density when scaled 
to Milky Way parameters; in addition, pd enters in the expression for Vh below to the 1/3 power. 
Thus, unless the fraction of collapsed baryons is unreasonably low, a value lower than unity does 
not lead to qualitatively different conclusions. The explicit expression for the blowout velocity, Ufe, 
has been obtained by Ferrara & Tolstoy (2000): 

where L is the mechanical luminosity of the parent SB. Note that v\, by definition (Ferrara &: 
Tolstoy 2000) is estimated at a height equal to 3H above the galactic plane. Obviously, if the 
blowout velocity has already decreased below the local sound speed at that point, the blowout 
is inhibited. By using eqs. ^, ^ and ^ it is easy to show that H^pd oc ^bC^, that is the product 
H^Pd is independent of mass and redshift; therefore Vh depends on these quantities only through 
L. Blowout occurs if Vb > Ve, or 



This implies that the mechanical luminosity for a SB to blowout must be larger than the critical 
value 



L, = {Q.27pGY''pdH' ( . (9) 
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The total mechanical luminosity of the corresponding SB in a primordial galaxy that virialized at 
redshift z can be written as (Ferrara 1998) 

Lt{z) = eouM^ = " '/' Mh, (10) 

where eo = 10^^ ergs is the energy of a SN explosion and is the star formation rate; we assume 
a Salpeter IMF, for which one SN is produced for each 56 Mq = of stars formed. The baryon 
density parameter is Qb = 0.05Of,^5, of which a fraction ~ O.OS/b^g (Ciardi et al. 2000a) is able 
to cool and become available to form stars. The free-fall time is tff = (47rG/?/j)~^/'^; = 3% is 
the star formation efficiency, which is fixed by matching to the cosmic star formation history as 
described in §4.2 below. 

For Population III objects, which as we have argued are likely to have been of low mass, 
we make the simplifying assumption that star formation is confined to a relatively small region, 
and that all the exploding SNae associated with the initial star formation episode result in the 
formation of a single superbubble. In this case a comparison between eq. ^ and eq. 10 shows that 



Lt{z)/Lc does not depend on mass and redshift, and the condition Lt{z) > Lc is satisfied when 

>0.1c^, (11) 

r 

which for typical parameters is equivalent to 4.4 x 10^^ > O.lCg, and is therefore valid for galaxies 
with Cs < 70 km/s. Thus low mass galaxies with coherent star formation (i.e. with all SNae 
driving a single SB) are very likely blowing out. 

For more massive and larger galaxies we have to consider the more realistic situation where 
the SNae are more widely distributed within the disk, and occur in OB associations with different 
values of Lqb (or, equivalently, with a different number N of SNae). In nearby galaxies it is found 
that the luminosity function of OB associations is well approximated by a power-law 

^(N) = ^= AN-^ (12) 

with /5 ~ 2 (McKee &: Williams 1997; Oey & Clarke 1998). Here Mqb is the number of associations 
containing N OB stars; normalization of 4'{N) to unity requires A = 1. Thus the probability for a 
cluster of OB stars to host N SNae is oc N~'^, where = LoBtos/^o, and toB = 40 Myr is the 
time at which the lowest mass (~ 8Mq) SN progenitors expire. The total mechanical luminosity, 
which must be equal to Lt{z) in eq. |lO|, is then found to be 

Nm 

Lt{z)= I LoB{N)ct>dN, (13) 

where Nm = 1 {Nm) is the minimum (maximum) possible number of SNae in a cluster. This gives 

Lt{z) = const. ^In^. (14) 

toB Nm 
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The contribution to the total luminosity from clusters powerful enough to lead to blowout is 

Lsiz, > Lc) = const. -^In^^, (15) 

toB J-^c 

where Nc is the number of SNae in a cluster with mechanical luminosity equal to Lc, i.e. 

iVe = (16) 

Thus, the fraction of the mechanical energy which can be blown out is 

Clearly, Nm (and therefore 5b) is an intrinsically stochastic number. To determine its dependence 
on the total number of supernovae Nt = Lt{z)toB / produced by a galaxy during the lifetime of 
an OB association, we have used a Monte Carlo procedure applied to the distribution function in 
eq. 12. The results for Nm and 5b as a function of Nt are shown in Fig. |l|; we recall that the star 



formation rate in the galaxy is related (through eq. 10) to Nt by M^, ~ 5 x lO^^A't Mq yr^^. 
As can be seen from Fig. |l|, for low values of Nt the quantity Nj^i is larger than N^ implying that 
in every galaxy at least some SBs are able to blow out. However, near Nt = 10^ Nm flattens and 
eventually becomes equal to N^ at A''^ ~ 45 000. Above this limit (corresponding to a galaxy with 
ili* K, 0.35 Mq yr"^ or Mh ~ 10^^(1 + zY^/'^Mq) blowout is inhibited. The fraction 5b can be 
seen from Fig. || to be a decreasing function of Nt, an approximate analytical form is 

5B{Nt) = 1 for Nt < 100 (18) 
5BiNt) = a + b ln{Nt'^) for Nt > 100, 

with a = 1.76, b = 0.165. Clearly, in small galaxies even the smallest associations are capable of 
producing blowout, so that the issue of coherence discussed above is irrelevant. 

Thus, if we assume the most of the baryonic matter has collapsed into the disk, the main 
metal polluters of the IGM at the present epoch are galaxies with visible mass lower than 

= QbMh ;S 5 X 1O^°M0. This critical mass decreases with increasing redshift as (1 + z) ' , 
indicating that the initial metal enrichment of the universe must have been produced by even 
smaller galaxies. For example, at z = 20, the expected value of M^^ is only 5 x 10^ Mq. 



2.2. Confinement of Blowout-driven Outflows 

Blowout will drive an outflow which will eventually be confined by the IGM pressure. What 
is the characteristic length. Re, at which such pressure equilibrium is achieved? By requiring that 
the outflow ram pressure, Pwv"^, is equal to the IGM pressure, pi we obtain 



- 8 - 



where we have used mass conservation and assumed that the flow is approximately spherical. Once 
recombination is complete at redshift z ~ 200, and before the IGM is reheated by the energy input 
from the first galaxies, the IGM pressure will evolve almost adiabatically, pi = + -2)/200)]^, 
where p-k/ks ~ 1500/i^ cm~^ K is the value of at 1 + z = 200. However, almost unavoidably, 
even the earliest SN-driven bubbles will expand in an IGM which has been pre-ionized by the 
same massive stars which later exploded as SNae. Thus, it seems more appropriate to calculate 
the IGM pressure in the surroundings of the galaxy as pi = {R/fi)pi{z)T, with T ~ 2 x 10^ K 
as a result of photoionization heating. This is a reasonable hypothesis since it can be shown 
Ciardi et al. 2000a) that the ionization spheres are much larger than the metal-enriched bubbles 
considered here. The mass loss rate is typically a fraction ^ ~ 10% of the star formation rate 
M^: = Le/eQV (Ferrara &: Tolstoy 2000) where = SsLt is the effective mechanical luminosity, 
that is the fraction available for blowout; we can also assume that Vw — Ve- With these hypotheses, 
and using the relations derived above for v^, one obtains 

1 /2 

R,(M,, z)= 0.9 (^^] cm, (20) 

\ -^4 Pi J 

where A4 = A/0.04. The behavior of Re with is shown in Fig. |2| for different redshifts. At 
lower redshifts, the bubbles are larger because the pressure of the confining IGM is lower; also 
evident from the figure is the increase of the critical mass for blowout with decreasing z. 

The major conclusion that can be drawn from Fig. ^ is that supernova-driven outflows are 
quite inadequate for dispersing heavy elements far from their production sites and cannot account, 
by themselves, for the ubiquitous presence of metals in the IGM at z ~ 3. This can be readily 
realized when we consider that in standard Cold Dark Matter models, taken here as representative 
of a larger class of hierarchical models of structure formation, the typical (physical) separation 
between objects with mass ~ IO^^Mq varies between 2 — 0.2 Mpc in the interval < z < 8 
(Ciardi et al. 2000a). This scale is much larger than the values of Re in Fig. |2[ so that we would 
expect regions outside the metal enriched bubbles to retain their primordial composition (this 
argument is reconsidered more rigorously in §4 below). Clearly some other mechanism must be 
operating to remove the metals from the surroundings of galaxies and mix them with the more 
generally distributed IGM. 

3. Predictions for CDM Models 

To make further progress, we need to fix a specific cosmological model. As an example, we 
consider the so-called Standard Cold Dark Matter (SCDM), with ^Im = 1,^A = 0, and h = 0.5; 
the power spectrum \6k\'^ is taken from Efstathiou et al. (1992), normalized to the present-day 
abundance of rich clusters {ag = 0.6; Eke, Cole, &: Frenk 1996). To calculate the number density 
of dark matter halos as a function of redshift we use the Press & Schechter (1974, hereafter 
PS) formalism; this technique is widely used in semi-analytical models of galaxy formation and 
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gravitational lensing (White & Frenk 1991, Kauffman 1995, Ciardi & Ferrara 1997, Baugh et 
al. 1998, Guiderdoni et al. 1998, Marri & Ferrara 1998) and it has been shown to be in surprisingly 
good agreement with the results from N-body numerical simulations. Given a power spectrum 
one can write the Gaussian variance of the fluctuations on the mass scale M: 

'r'M = j ^,W\k,R)\8,\\ (21) 

where M = (4/3)7r/>i?^, p is the matter density, and 

W = — ^ [sm{kR) - (kR) cos{kR)] ; (22) 

is a top-hat filter function. From the results of the nonlinear theory of gravitational collapse, 
stating that a spherical perturbation with overdensity 5c = Sp/p > 1.69 with respect to the 
background matter collapses to form a bound object, through the PS formalism we can derive the 
normalized fraction of collapsed objects per unit mass at a given redshift: 



f{M, z) = jlMi±i) e-^?(i+^)V2.I, (_<^\ . (23) 
\ TT alj \ dM J 

Then the comoving number density of dark matter halos per unit mass is 

n,{M,,z) = ^^^f{M,,z). (24) 

We can now ask what is the fraction of the IGM polluted by metals at different redshifts. To this 
end we calculate an IGM porosity parameter, Q{z), defined by 

dQ{z) = Mh\^\R\Mh,z)dz, (25) 

in two different cases. First, we consider the case in which metals are only injected in the IGM 
by superbubbles and therefore R = Re in the previous equation. The resulting porosity evolution 
(complete overlap of the bubbles occurs for Q = 0.16, Smith 1976) is shown in Fig. ^. Again we 
see that blowout by itself would lead only to a negligible dispersal of metals {Q < 10^^), with 
most of the IGM maintaining its primordial composition to the present day. 

If we are to explain the relatively ubiquitous presence of metals in the 'true' IGM, as deduced 
from at least some observations, we are then forced to assume that some additional physical 
mechanism, capable of efficiently transporting metals away from their production sites, must be 
at work. The nature of such a process can only be matter of speculation at present, because 
neither the numerical simulations nor the observations have yet reached the required levels of 
sophistication or sensitivity to address this question properly. Among the options which are 
worth exploring are galaxy collisions, diffusive processes and the peculiar motions of galaxies. We 
consider the relative importance of these different possibilities in future work. 
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Here we take a strictly phenomenological approach which nevertheless has considerable 
predicting power. We introduce a diffusive radius, Rd, defined as the mean interdistance between 
the galaxies responsible for providing the predominant contribution to the metal enrichment of 
the IGM (see Fig. g), i.e. = 0.2 Mpc ( comoving) for > 2 x Wi^Mq. As the largest Doppler 
parameters measured in Lya clouds are of order ~ 50 km s~^, the time required for a pollution 
front to travel such a distance is shorter than the Hubble time only at redshifts z <1. Thus, our 
assumption is equivalent to fixing at z ^ 1 the epoch at which metals become homogeneously 
distributed. As we will see shortly, this simple hypothesis leads to a number of consequences 
which are in accord with available data. If future studies of possible mechanisms for the diffusion 
of metals are able to determine directly the value of Rd and its evolution with time, it will be 
relatively easy to include their results in the general framework of this paper and explore any 
differences with the conclusions presented here. 

As can be seen by Fig. ^, when the condition R = max{R(,, Rd) is introduced, mixing is 
improved dramatically — as expected, and metal-enriched bubbles indeed overlap at z ~ 1. At 
later epochs essentially all of the IGM is polluted with heavy elements produced in galaxies and 
subsequently redistributed by the combined effects of blowout and diffusion. 



3.1. Metallicity of Polluted Regions 

We can now calculate the average metallicity of the gas inside the diffusive spheres. In 
order to do so we need to know fiz, the mass of metals produced by the typical supernova. 
Nucleosynthesis calculations in Type II SNae by Tsujimoto et al. (1995) predict /i^ « 2 — 3Mq 
for a Salpeter IMF, depending on the upper mass limit above which a black hole is formed. Here 
we adopt nz = 2.58Mq, also consistent with a matching to the cosmic star formation history 
(§4.2). By design, mixing is efficient and the metal distribution is therefore homogeneous within 
the radius Rd- The total number of SNae per galaxy is 

MsN = ^M,. (26) 

Then the mass of metals ejected by a galaxy inside a halo of mass is 

Me = t^^^M,6s, (27) 



with 5b given by eq. |18|. The metal density in a diffusive sphere surrounding a given halo is 
pz = 3M(,/4:TtR^, and its average metallicity is pz/pi{z), where Pi{z) is the mean density of IGM 
at redshift z. This assumption restricts our analysis to the case of the 'true' IGM, characterized 
by overdensities close to unity. Extrapolation of the results to hydrogen column densities 
logNffj ^ 14 is only very qualitative and deserves more study. The dependence of metallicity 
in the diffusive spheres on halo mass and redshift is plotted in Fig. ^. At any given redshift, Z 
increases as a function of mass (because of the increasing metal production) up to an abrupt cutoff 
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when blowout becomes inhibited, as described by the behavior of Sb in Fig. ||. As the spheres 
grow with time, the metals are distributed over larger volumes and this is reflected by the shift to 
lower metallicities with decreasing z in Fig. ^ Note that for the same reason pockets of very high 
metal content {Z ~ l-^©) are expected at large z, although the size of these regions is very small. 

3.2. IGM Metallicity Evolution 

Once the cosmological model has been fixed, we need to determine the values of the the star 
formation efficiency, r^, in order to calculate the metallicity of the IGM as a function of time. We 
choose this value by comparison with observations, as follows. We calculate the evolution of il^, 
the density parameter of stars formed in the universe (see Fig. ^), with the requirement that at 
z = it matches the estimate by Fukugita, Hogan &; Peebles (1998). These authors conclude that 
il^(O) in spheroids, disks and irregulars is ~ O.OO49/15Q. A similar value is found by integrating 
current estimates of the star formation rate density in the universe as a function of redshift, as 
deduced from deep imaging surveys (Pettini 1999). With this normalization (giving « 3%), 
we can then derive the corresponding evolution of the metals produced by stars and returned to 
the ISM (Fig. |5|), which is directly proportional Q^,. 

However, not all the metals can escape from the galaxy where they have been produced. The 
cosmic ejection fraction, fej{z), (i.e. the ejection fraction averaged over the entire population of 
halos) is very close to unity at high redshift where predominantly small galaxies are present, but 
it steadily decreases to about 50% at z = 0, as the number of more massive galaxies able to retain 
their metals increases. For this reason, the curve describing the density parameter of ejected 
metals = fej{z)Qz{z) in Figure 5 progressively deviates from that for Qziz) with decreasing 
redshift. Stated differently, today about 50% of the metals produced should reside in the IGM. 
If the metals were homogeneously mixed with the baryons in the universe at any redshift the 
average IGM metallicity (top curve in Fig. I) would be {Z) = /Qb - 1/25Zq at z = 3 and 
{Z) ~ O.IZq at z = 0. Note that the average metallicity of today's galaxies would be higher by a 
factor r2fe/(0^ + Qg) (the baryon density divided by the sum of star and gas density in galaxies); we 
therefore expect that for luminous matter today {Z) Zq, in accord with observational estimates 
(Edmunds & Philhps 1997). 

Our main results are displayed in Figure 6 and 7, which show in the spread in metallicity as a 
function of redshift, compared with the average {Z) of the IGM from Figure 5. The two figures 
correspond to different metallicities of the outflowing hot gas in the superbubbles, assumed to be 
Zq in Fig. 6 and 8Zq in Fig 7. In the figures the metallicity distribution at each redshift considered 
is shown in the vertical direction with the density of symbols approximately proportional to the 
amplitude of the distribution. The distributions are relatively flat, but at redshifts z ^ 5 they also 
show a double-horned profile with points accumulating at the highest and lowest metallicities. We 
illustrate this effect in Figure 8, where we have reproduced the metallicity histograms at 2; = 
and z = 1 from Figure 6. The more massive galaxies are responsible for the high metallicity peaks 
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in these distributions, caused by the cutoff in the ejection fraction described by the parameter Sb 
(Figure 1). The low metalhcity maximum is due to the increasing number of low mass halos. 

There are several interesting features of Figures 6 and 7 which we now discuss briefly. First, 
at z > 6 the average IGM metallicity is not within the predicted ranges of metallicities, being 
lower than the lowest values in the distributions. What we are seeing is a highly inhomogeneous 
distribution of metals which are still clumped in small regions around the galaxies which produced 
them. As mixing proceeds, the average metallicity increasingly becomes a better description of 
the true mean metallicity, and approaches the mean of the distribution. At z ~ 1 the postulated 
diffusion leads to all the volume in the universe being polluted with heavy elements to some extent 
(Figure 3). The merging of the metal-enriched spheres produced by different galaxies has the 
effect of erasing regions of low metallicity: the minimum metallicity corresponds to that in spheres 
marginally overlapping at a given redshift. This effect is reflected in the growth of the lower 
boundary of the distributions from (1 + z) ^5 in Figures 6 and 7, and results in the present-day 
IGM metallicity being confined within the narrow range 0.1 it 0.03Zq(1 a). 

The decreasing spread in metallicity with time is a direct result of the growth of {Z). At high 
z, where the overall level of enrichment is low, even relatively metal-poor spheres have a chance to 
stand out. However, as the average IGM metallicity increases, only a few massive objects are able 
to produce diffusive spheres with Z sufficiently high to be recognized as metallicity peaks (and 
still be able to eject their metals). The situation is similar to a "flooding effect", where as the 
water level rises fewer and fewer mountain peaks can be seen. 

At redshifts z > 1, when metal enrichment is highly inhomogeneous, a considerable fraction 
of the IGM is still of essentially primordial composition. Thus, at z ~ 3 for example, we expect 
that only some low column density Lya forest clouds will show associated metal lines, while the 
majority will not; the ratio between Lya clouds with and without metals depends on the covering 
factor of the diffusive spheres and grows with time. This has implications for the interpretation of 
the results of searches for metals in the forest. Given the low optical depths when log A^(II I) < 14, 
such searches are normally conducted in a statistical way, by considering together the data from 
many absorption lines. In such cases one obtains some gross average over all the absorbers 
which, given the dilution with the unpolluted IGM, will be systematically lower than the values 
of (Z) plotted in Figures 6 and 7. For this reason we also show in these figures the covering 
factor-weighted metallicity {Z P{z)), where P{z) describes the evolution of the covering factor 
down to z = 1 . 

In Figures 6 and 7 we have used different symbols to represent the contributions to the IGM 
metallicity from halos with virial temperatures above (hexagons) and below (triangles) 10^ K, to 
give a qualitative idea of the role of large and small objects in the enrichment process. It is seen 
that at high z small objects are controlling the process, whereas at lower redshift enrichment is 
largely regulated by more massive galaxies (up to the ejection limit set hy 6b, see Figure ^). If 
some inhibiting effect, such as photoheating by the UV background, affects preferentially low mass 
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galaxies, the lower metallicity bounds at high redshift will move up accordingly, shifting closer to 
the line separating hexagons from triangles. The flat upper boundary of the distribution at high 
z reflects the fact that the volumes involved are so small that the IGM baryon loading results in 
a negligible dilution of the metallicity of the ejecta. Since there is no firm measurement of this 
quantity at present (Heckman, private communication), we have considered two possibilities, Zq 
(our standard case) and 8 Zq in Figures 6 and 7 respectively. The IGM metallicity distribution is 
affected by this choice at high redshift, but the difference becomes much smaller at low z, where 
the IGM baryon loading regulates the dilution of the larger diffusive spheres. 

Finally, we also show in the two figures the average metallicities for the case in which 
metal-enriched spheres with size below Re{min) = 1 kpc have been excluded. In our models these 
are the absolute lower limits to the average IGM metallicity at a given redshift. R(,{min) has 
been calculated by imposing the condition that the collisional timescale between galaxies with a 
given impact parameter is shorter than the Hubble time at any redshift up to z = 10. Obviously, 
a sphere with size greater than Rg{min) can be strongly disturbed by tidal interactions following 
an encounter at greater impact parameter, and for this reason the curves obtained in this way are 
strictly lower limits. Nevertheless, as can be seen from Figure 6 and 7, imposing this condition 
reduces the metallicity by only about 40%. 

In addition to the SCDM model, we have also considered a CDM model with a cosmological 
constant (CDM+A) with Qm = 0.4, = 0.6, /i = 0.5). By normalizing the star formation rate 
following the same procedure as for the SCDM model we obtain a higher value of = 10.5%. 
Because of the normalization, the resulting metal distribution is qualitatively very similar to 
the one derived for the SCDM case. The only notable differences are a larger metallicity spread 
between 0.5 < z < 2 and a somewhat higher mean value at z = 0, (Z) = 0.15 it O.OSZq. 

4. Discussion 

The basic conclusion of this paper is that metal ejection driven by SN events fails, by more 
than one order of magnitude, to distribute the products of stellar nucleosynthesis over volumes 
large enough to pollute the whole IGM to the typical metallicity of Lya clouds, [C/H] ~ —2.5. We 
are therefore forced to conclude that some additional physical process must be at play, the nature 
of which remains to be determined. In our scheme transport of metals occurs in two sequential 
steps. First, SNae provide the initial kick necessary to eject heavy elements outside the potential 
well of a galaxy, but the ejecta are then confined by the IGM pressure to a relatively small bubble, 
of radius Re- We then postulate that a second process is responsible for the diffusion of metals 
on a typical scale R^, comparable to the mean distance between the galaxies which are the most 
efficient pollutants of the IGM. 

Our scenario is qualitatively different from that proposed by Gnedin (1998), who attributed 
the mixing to more violent and rarer galaxy mergers. In the model proposed here SNae are of 
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fundamental importance as they initiate the mixing process. However, galaxy collisions might well 
play a role in the subsequent phase during which metals are spread over larger scales, of order i?^. 
Other diffusive processes, such as thermal conduction and turbulent mixing layers at the interfaces 
between cosmological flows, as well as the peculiar motions of galaxies, may also be important in 
determining the spatial structure of the distribution of metals. A detailed study of such processes 
is beyond the scope of this paper and should, in any case, be based on cosmological simulations 
which are in progress. We will report of this extension of the work in a forthcoming paper. 

In the present study we have stressed that the IGM enrichment proceeds in a very 
inhomogeneous manner, with pockets of metal rich gas gradually increasing both in number and in 
size until they overlap at 2; ~ 1. The average metallicity of the IGM increases with time, a trend 
confirmed by the results in Barlow & Tytler (1998), who found an order of magnitude increase in 
the metallicity of the Lya forest between z = 2.5 and z = 0.5. Although our results are strictly 
applicable only to the 'true' IGM with overdensities close to unity, we nevertheless regard this as 
an encouraging performance of our model. 

The metallicity spread is predicted to decrease with the progress of time. At z ^ 1, when the 
entire volume of the universe has been exposed to metal pollution, the spread in metallicity is less 
than one order of magnitude, and at the present epoch [Z) = 0.1 ± O.OSZ©. Thus we predict that 
at 2: < 1 essentially all absorbers should have associated C IV absorption, irrespectively of their 
column density. The effect should be very pronounced. Not only are the voids polluted by the 
overlap of metal-enriched diffusive spheres, as discussed above, but the decreasing intensity and 
hardness of the ionizing background lead to a further increase in the fraction of C which is triply 
ionized, so that the ratio N{C IV)/7V(H I) increases for a fixed [C/H] (Ranch et al. 1997). It 
should be possible to test these predictions with forthcoming observations. UV-efHcient echelle 
spectrographs, now available on the VLT and soon on the Gemini South telescope, will allow 
sensitive searches for C IV doublets at significantly lower redshifts than probed so far, while STIS 
on EST and FUSE will map the Lya forest with the required spectral resolution at wavelengths 
below 3000 A, which are inaccessible from the ground. 

While our study has concentrated on the IGM, it also has important consequences for the 
metallicity of the intracluster medium, which is found to have a rather uniform value k. I/'^Zq 
with little, if any, evolution up to 2; 2± 0.3 (Fukazawa et al. 1998, Renzini 1999). This has been 
interpreted as evidence supporting the view that the enrichment process in clusters was already 
completed by that epoch. In general we expect two sources to contribute to the build up of the 
intracluster medium, infalling IGM and gas stripped from cluster galaxies. In our models these 
two components have significantly different composition. If the clusters formed at 2; ~ 1, we 
expect the IGM to have O.IZ© (roughly constant from z ~ 1 to the present), while the galaxies' 
ISM has approximately solar composition. The simplest estimate of the resulting metallicity is 
the geometric mean of these two values (as appropriate for hydrodynamical mixing problems, see 
Begelman & Fabian 1990), that is {Z)icm = (0.1 x IY^'^Zq = 0.32^©, as observed. While this may 
well be a fortunate coincidence, it is also true that this conclusion lends qualitative support to the 
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model and the assumptions made. 

We make a final point concerning the role of intergalactic dust. The gas phase abundances 
we have derived do not take into account the possibility that some fraction of heavy elements 
may be locked up into dust grains. If this were the case, clearly the distributions of metallicities 
in Figure 6 and 7 would shift towards lower values of Z . In addition, at redshifts were the 
enrichment is still inhomogeneous, dust associated with regions of high metal concentration may 
reprocess UV/optical photons into IR radiation and give rise small scale anisotropics in the Cosmic 
Microwave Background which may be detectable (Ferrara et al. 1999). Having said this, however, 
we consider it highly speculative whether dust can survive at all in the hostile environments 
associated with outflowing superbubbles, where the temperatures are high and the gas has been 
shocked by SN explosions. 

5. Summary 

In this paper we have investigated the evolution of the metallicity of the intergalactic 
medium with particular emphasis on its spatial distribution. We have derived the conditions 
under which supernova-driven ejection of metals from galaxies can occur. A strong conclusion of 
our calculations is that if SNae were the only source of kinetic energy for the metals, a highly 
inhomogeneous distribution would result at any redshift. Under these conditions we would expect 
most of the volume of the universe to remain at near-primordial composition, with a metallicity 
Z ~ W^'^Zq, in contrast with the observational results discussed in the Introduction. Thus, an 
additional (but yet unknown) physical mechanism must be invoked to mix the metals on scales 
comparable to the mean distance between the galaxies which are the most efficient pollutants. 
From this simple hypothesis we have derived a number of testable predictions for the evolution of 
the IGM metalhcity. 

Quantitatively, we find that: 

1. Metal ejection, or blowout, is inhibited in galaxies with total mass above 
M/j 10^^(1 -I- z)~^/'^Mq due to the combined effects of their larger gravitational field 
and less coherent SN energy deposition. 

2. The fraction of metals ejected over the star formation history of the universe is about 50% 
at z = 0. We expect that at the present epoch approximately half of the metals are to be found in 
the IGM and the average metallicity of luminous matter to be approximately solar. 

3. If the ejected metals were homogeneously mixed with the baryons in the universe, the 
average IGM metallicity would be {Z) = rj^/r^;, ~ 1/25^© at z = 3. However, due to the 
spatial inhomogeneity, {Z) is actually lower than the mean of the distribution in the diffusive 
metal-enriched spheres. 

4. Metals become homogeneously distributed in the IGM at z ^ 1, when the metal-enriched 
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zones overlap, and the spread of the distribution is reduced. We calculate that at z = the IGM 
metallicity is in the range Z ~ 0.1 it O.OSZq. 

5. The uniform metal abundance of intracluster gas from z ~ 0.3 to the present is naturally 
explained by a mixture of IGM infalling into the cluster and gas stripped from cluster members, 
with element abundances for both components as predicted by our models. 

We should like to thank E. Corbelli, A. Meiksin and B. Nath for useful discussions. YS 
acknowledges support from the OAArcetri. 
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Fig. 1. — Upper panel: dependence of the maximum number of SNae in a cluster as a function of the total number 
of SNae, Nt, produced by a given galaxy in an OB association lifetime as obtained from Monte Carlo simulations; 
also shown is the critical number of SNae, A^c, for blowout. Lower panel: Energy fraction available for blowout, &b 
as a function of the total number of SNae produced by a given galaxy in an OB association lifetime as obtained from 
Monte Carlo simulations; also shown is an analytical fit to the numerical data. 
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Fig. 2. — Radius of the metal bubble (in physical units) as a function of the halo mass for different virialization 
redshifts z — 0,1,2,3,4,5 from the uppermost to the lowermost curve, respectively. The endpoint of each curve 
denotes the critical mass for blowout. 



- 21 - 




Fig. 3. — Porosity factor of the metal bubbles considering transport by ejection driven by superbubbles (lower 
curve) and transport due ejection and diffusion processes acting in combination (upper curve). Overlapping of the 
various bubbles is seen to occur at z « 1, where Q{z) is equal to the critical value 0.16. 
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Fig. 4. — Metallicity of diffusive spfieres as a function of parent halo mass for redshifts z — 0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4 
from the bottom to the top curves, respectively. The curves shown are for the case where the initial metallicity of 
the hot gas in the superbubbles is solar. 



- 23 - 




- 24 - 



10^ w 



0.1 



ISI 



10 



-2 



10 



-3 



10-* r 



T r 



T 1 1 — I =3 



iZ(R,>l kpc)/<Z> 




J I I I I I L 



10^ 



1+z 



Fig. 6. — Redshift evolution of the IGM metallicity distribution in the SCDM model. Also shown are the average 
IGM metallicity, (Z), the contribution from objects with large metal-enriched spheres, Z{Re > 1 kpc), the ratio of 
the two, and the covering factor-weighted metallicity, {ZP{z)). The initial value of the hot gas metallicity is Zq 
(standard case). 
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Fig. 7. — Same as Fig. 0, but for initial value of the hot gas metallicity 8 Zq. 
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Fig. 8. — IGM metallicity distribution function for the SCDM model at redshift z — 1 (obliquely dashed) and z = 
(solid) . Both histograms are normalized to unity; the initial value of the hot gas metallicity is Zq . 



